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Abstract 

We investigate the evolution of a system composed of A'^ non-interacting point particles of 
mass m in a container divided into two chambers by a movable adiabatic piston of mass M S> m. 
Using a two-time-scale perturbation approach in terms of the small parameter a = 2m/ (M -|- m), 
we show that the evolution towards thermal equilibrium proceeds in two stages. The first stage 
is a fast, deterministic, adiabatic relaxation towards mechanical equilibrium. The second stage, 
which takes place at times 0{M), is a slow fluctuation-driven, diathermic relaxation towards 
thermal equilibrium. A very simple equation is derived which shows that in the second stage, the 
position of the piston is given by Aai (t) = L[l/2 — ^{cet)] where the function ^ is independent 
of M. Numerical simulations support the assumptions underlying our analytical derivations and 
illustrate the large mass range in which the picture holds. 
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1 Introduction 

The "adiabatic" piston problem is a well-known controversial example of thermodynamics [1]. An 
isolated cylinder contains two identical fluids which initially are in different equilibrium states and 
which are separated by an adiabatic fixed piston. The whole system remains therefore in equilibrium. 
The problem is then to predict the final state to which the system will evolve when the constraint 
fixing the piston is released. Although it is a very old problem since it was discussed to measure 
experimentally the ratio Cp/c^ of gases already before 1940 [2] [3], it still remains a controversial 
question since it shows that the two laws of equilibrium thermodynamics are not sufficient to predict 
the final state. 

Recently, especially after the work on the second law of thermodynamics by E. Lieb and J. Yn- 
vagson [4] and the talk of E. Lieb at the StatPhys 20 meeting in 1998 [5], this problem has attracted 
renewed interest. It was first realized that it is a standard example where one is forced to apply 
non-equilibrium thermodynamics since the final state may depend on the values of the friction coeffi- 
cients [6] . Then a very simple model was considered to investigate the evolution from a microscopical 
point of view [7] . This system consists of N non-interacting particles in a cylinder of length L and 
cross-section A. It is divided into two compartments containing respectively N~ and particles 
(A^ = N~ + N^) by an adiabatic (i.e. no internal degrees of freedom) piston of mass M. The dynamics 
is defined by the condition that the piston is constrained to move without friction along the x-axis 
and the particles make purely elastic collisions on the boundaries and on the piston. Without loss 
of generality, we can assume that all the particles have velocities parallel to the x-axis and thus we 
are led formally to a one-dimensional problem (except for normalization). Therefore, given that the 
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velocities of a particle and the piston are v and V before they collide, then after the collision the 
velocities will be v' and V, with: 



v' = 2V -v + a(v-V) V' = V + a(v - V) where a = — (1) 

For physical situations where m M , this model was investigated in [8] using Boltzmann equation 
and a perturbation expansion of the velocity distribution function of the piston in powers 

of e = y/m/M. For an infinite cylinder (L — oo) and equal pressures on both sides of the piston 
{p~ = p^). it was shown that the stationary solution of the Boltzmann equation gives a constant 
velocity for the piston Vgt = (1/4M) y/2TrkBm {^/T+ — ^/T^)+0{m/M) towards the high temperature 
domain. It was thus concluded that stochastic motion together with space asymmetry (temperature 
difference) implies a macroscopic motion. Using qualitative arguments and numerical simulations, the 
case of a finite cylinder was investigated in [9]. It was thus realized that the evolution takes place 
in two or three stages. In a first stage, the evolution appears to be deterministic and adiabatic; this 
first stage proceeds until mechanical equilibrium is reached = p"*" but T~ ^ In the second 

stage which takes place on much longer time scale, the simulations showed that the motion of the 
piston is stochastic and proceeds with exchange of heat through the piston until thermal equilibrium 
is reached where T~ = r+, i.e. the piston which was adiabatic when fixed and during the first stage 
becomes heat-conducting under fluctuations. This explains in particular the results mentioned above 
for the infinite cylinder with p~ = p+. In the third stage, under the stochastic motion of the piston, 
the velocity distribution functions of the fluids tend to Maxwellian distributions. 

Recently we have studied the adiabatic piston problem in the thermodynamic limit for the piston, 
i.e. by considering the limit where L is fixed but the area A of the cylinder, the mass M of the piston 
and the number N"^ of fluid particles tend to infinity with M/A and = mN^/M flxed [10]. Starting 
from Liouville equation, it was shown that in this thermodynamic limit, the motion of the piston is 
adiabatic and deterministic, i.e. (l^")f = {V)^. Introducing at this point simplifying assumptions 
(see Assumptions 2 and 3 below), we obtained a system of autonomous equations from which we 
concluded that the system evolves towards a state of mechanical equilibrium where the pressures 
are equal but the temperatures different. Furthermore, numerical simulations were conducted which 
showed that the motion depends strongly on for < 1 but tends to be independent of for 
> 10. For R^ < 1, the evolution is very weakly damped and the period of oscillations, which 
depends on i?^, coincide with the period computed with our equations, and with the period obtained 
from thermodynamics assuming adiabatic oscillations [2] [3]. On the other hand for R^ > 10, the 
motion is strongly damped and is independent of R^. An equation was thus conjectured to describe 
the evolution for R^ large enough. It should be noticed that these two types of regimes, i.e. weak vs 
strong damping, have been observed experimentally [3], as well as in the previous numerical simulations 
for the simple piston where the evolution is necessarily adiabatic [11]. Let us also stress that in this 
thermodynamic limit, we could deduce from Liouville equation the factorization property of the joint 
velocity distribution for the piston and a particle at the piston surface. 

In the present work we investigate the motion of the piston with M finite but M ^ m. Introducing 
now the factorization property as an assumption (sec Assiunption 1 below) which we assume to be 
valid to first order in a, and using a two-time-scale perturbation approach, we show that the evolution 
of the piston proceeds in two stages with totally different properties and time scales. In the first stage, 
characterized by a time scale ti — Nm/ E^) (average time for a particulc to collide twice with the 
piston) where Eq is the initial energy, the evolution is adiabatic and deterministic, with properties 
similar to those obtained in the thermodynamic limit M = oo (i.e. a = 0). In particular, the evolution 
can be either weakly or strongly damped and for M sufficiently large the evolution is independent 
of M (for fixed i?^). This evolution proceeds until "mechanical" equilibrium is reached where the 
pressures become equal and it lasts for a time length of order 0{M/m). 

In the second stage, characterized by a time scale t2 ~ Mti /m, the evolution is strongly dependent 
on M. In fact, we shall derive the following scaling property: 
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Main result: introducing the scaled time variable t = at, the evolution of the piston of mass M ^ m 
is given by: 

'1 



XM{r) = L(--ar)] (2) 



where 



dT 



1 


1 
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L ^ 







TV / 



(3) 



Moreover from the knowledge of ^(r), one obtains the temperatures and pressures on both sides of 
the piston, as well as the temperature of the piston, to lowest significant order in a. In this second 
stage, the fluctuation-driven evolution is stochastic and proceeds with heat transfer across the piston 
towards a state of "thermal" equilibrium where T~ = T+. Investigation to higher orders in a would 
be necessary to conclude that the fluid velocity distributions will ultimately tend to Maxwellian 
distributions. 

The same microscopical model was analyzed by Lebowitz, Piasecki and Sinai using a different 
limiting procedure [12]. In their work, they have considered the container to be a cube of size L and 
they have taken the limit L — » oo with A''" = ~ and M ^ L^. Using heuristical arguments, 
they derived autonomous coupled equations describing the motion of the piston and the fluid for large 
L. More recently, exact results were obtained by Chernov, Lebowitz and Sinai using this limiting 
procedure L — » oo [13]. They were able to prove rigorously that for a time interval sufficiently 
short so that only the first and second recollisions of a particle on the piston can occur, the random 
functions describing the position and the velocity of the piston, expressed in terms of scaled variables 
(r = t/M, Y = X/M), converge in probability to some deterministic functions. 

To conclude this introduction, we should mention that the approach to "thermal" equilibrium was 
investigated in [14] using an expansion of the master equation for the piston velocity distribution 
function in powers of ^/a, as well as numerical simulations, for a gas of hard rods. 

Another interesting discussion was presented in [15], where the power spectrum and the time 
correlation function of the piston were calculated at equilibrium for a two-dimensional system of hard 
disks by molecular dynamics simulations. Prom this analysis, one can also deduce a two-time-scale 
relaxation towards equilibrium. The problem of approach to equilibrium and the question of heat 
transfer for this 2-dimensional system of hard disks was analyzed in [16] and [17] using molecular 
dynamics simulations and a linear Langevin equation. 

2 Equations for the moments (V^) of the piston velocity 

In the spirit of BBGKY hierarchy, we shall characterize the velocity distribution ^{V,t) through its 
moments: 



m = {vm 



/oo 
V^V,t)dV (4) 
-CSO 



and 



/oo 
V"-^V,t)dV (5) 
-OO 

Introducing the notation: 



/OO rV 
^ - ^)'p;„r-/(«' V' t)dv - J_Jv - V)''pt^^f{v, V, t)dv 



(6) 



where A; > and /oJu^j(w, V,t) is the two-point correlation function for one molecule on the left (— ), 
resp. on the right (+), and the piston, then by integration of Liouville equation over all variables 
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except the piston velocity V, we have obtained the foUowing evolution equation for $ [10]; 



k=0 



where 7 = Aa = 2mA/ {M + m) and with dt^{V,t) defined in the framework of generalized functions. 
Hence, by definition, integration over V commute with dt and with the infinite summation in any 
average with respect to 

In particular, we have: 

1 _ (-i)^a^+^ , 



Under the hypothesis that V — > P^r/(^'^'*) decreases faster than any power of V, uniformly in v, 
we may integrate by parts, which yields: 



All terms with k > n, appearing as pure derivatives with respect to V, vanish after integration on the 
whole y-axis, due to the fast decrease ofV-^ pfur/(^' infinity. Hence the equation simplifies 

into: 



where depends functionally on Pg^^f{v, V,t). 



3 Factorization property 

In [10], we have shown that in the limit a = 0, the two-point correlation functions pf^rfi''^' ^) ^^^^ 
the factorization property. Although we can not prove it, we expect that in the limit a — > which 
we now consider, and for initial conditions such that the evolution is smooth, this property will still 
be valid to first order 0{a) where the leading behavior is ^{V,t) = 6{V — V(t)). Therefore in the 
following perturbation approach, we shall consider that the following assumption holds to first order 
in a: 

Assumption 1 (factorization property): the two-point correlation functions have the following factor- 
ization property at first order in a: 

= t) *(^' (12) 

Under this factorization property, we have = where: 

/oo nV 
^ i.v-Vfp-^^f{v,t)dv-j {v-Vfpt^^f{v,t)dv 

(13) 



4 



and the evolution of the moments (y") of the piston velocity satisfy the equation: 

Similarly, under the factorization assumption, the distributions p^{x, v, t) of the fluid molecules satisfy 
the Boltzmann equations with boundaries [10]: 

{dt +vdx)p~{x,v,t) = S{x + L)vp~{-L,v,t) 

+6{x - X{t))[V{t) - vMV{t) - v)p-{X{t),v',t)+e{v - V{t))p-{X{t),v,t)] (15) 

{dt + vdx)p'^{x,v,t) = -5{x - L)vp+{L,v,t) 

-S{x - XmV{t) - v][0{v - V{t))p+{X{t), v', t) + e(V(t) - v)p+{X{t), V, t)] (16) 

where v' = 2V — v + a{v — V) is the velocity of the molecules after their collision onto the piston. 
Note that p^{X{t),v,t) is what we denote pfurfi''^^^)- 

It is straighforward to check that: 

-^iV,pLf) = -kFt,{V,pLf) (if k > 1) (17) 

which leads to define for fc < by: 



(^)' ' fi..,(V.t) {r>l) (18) 



We define the densities p^^j (not to be confused with the distributions temperatures 
^sur/ ^^'^ pressures pf^^f surface of the piston, on each side: 



Pjurf = 2/ iv,t)dv = 2F,-{V = 0,p;^^^) (19) 

Jo 

Pturf = f pturA^,t)dv = 2F+{V = 0,pj^,^.) (20) 

J — oo 

Pfurf = '^rnF^iy = ^)=ptrfkBTf^rf (21) 

It is obvious from the definition that F^ {V) and — F2{y) are decreasing functions of V ^ which allows 
to express F^iV) and F^iy) in the following form: 

M + m\ ,± 



2mF^{V) = ± [—^ j A±(y)y (22) 

and A^(y) have the physical meaning of friction coefficients. Denoting X{V) = X~{y) + \^{V) yields: 



(24) 



From Eq. (17), we have also 2mF2{V) = p^ — 'imF^{0)V + 0{V'^), so that, denoting simply 
A± = A±(y = 0): 
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4 Irreducible moments 



We denote Fn' the r-th derivative of F„ (with respect to V); from Eq. (17), these are functional 



(r) 



of P^y^rf{v,t), proportional to Fn-r{V, p^^^^). We thus have: 
where: 



r=0 



(F±)W = (-l)-- 



(n-r)! 



for r < n 



(F±)M = ±(-irn! 



d_ 
dV 



r—n—1 



for r > n + 1 



Plugging the expansion Eq. (26) in Eq. (14) yields for n = 1: 



1 dV 
^~dt 

where: 

p(0) _ p. c.(l) _ OP. f(2) 



(26) 

(27) 
(28) 

(29) 



i^r=^2 F2^^^=-2Fi i^2'''=2i^o i^r^=2i^_i = -2[p,-,/F,t)-p+,/y,t)] (30) 
and Ar axe the irreducible moments: 



A, ^ r {V-V{t)r^{V,t)dV = ^(-1)^ ^'W'^ {V-'') 



with 



/oo 
y$(y,t)dy 
-oo 

the average velocity of the piston. Let us note that 



Ao = 1 
Ai = 

A2 = (y2)-y^(i) 



(31) 



(32) 



(33) 
(34) 
(35) 



We can invert Eq. (31) to express the moments (V^) as functions of the irreducible moments (cumulant 
expansion): 



si 



(36) 



Finally replacing the function Fk_^_2{V, pf^^^) by its expansion in powers of {V — V{t)) around the 
average velocity V{t) and using Eq. (36) lead to the following evolution equations for the irreducible 
moments (s > 2): 



1 dAs 



7 dt 



= -sAs_i 



r>2 



(r) 



2s 



(n+l)!''^+"'^i 



k=0 



{k + 2y.{s - 2 - ky. 



-2-fe+n-f3+fe 



n>0 



(37) 
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where we recall that Fk are functions of V and functional of p^{.,t). 
For s = 2, it comes: 

IdAo . A, 



4E7^^i'"' + -E>"^i"^ (38) 



7 rft 4-^, (n + 1)! ' ^ n! 

' n>0 ^ n>0 



Using the fact that Ai = 0, F^'^'^ = -5^^" and F^''^ = -3F^'' if r > 1, we may rewrite this 
equation: 

Y^^{2r-3a)F^^-'^ + aF, (39) 



1 dA2 ^ ^ ^O. _ f('-I) 

' r>2 



5 Evolution at first order in a 

A qualitative analysis of Eq. (37) shows that: 

As ~ al^l (40) 

where [{s + l)/2] is the integral part of (s + l)/2. Since we want to restrict our study to first order in 
a, assuming that Eq. (40) is valid, wc only have to consider the evolution of V and A2, Eqs. (29) and 
(39), which we supplement with the equations for the energies (E^) of the gas in the left and right 
compartments, restricted similarly to first order in a. We thus obtain to first order in a the following 
set of coupled deterministic equations: 



IdV 

\dA2 

7 dt 



F2 + A2F0 (41) 
-4A2J^i + aFs (42) 



= -My[F2-+A2Fo-] + f [4A2i^r-"^^3-] (43) 
1^ = MV[F} + /^.F+]-^[^l^,F+-aF+] (44) 



and we recall that all the functions F are functions of V{t) and functionals of pfuj.f{v,t). We should 
note that the set of Eqs. (41)-(44) have a constant of motion, i.e. 

{E-) + {E+) + ^ M{V^ + A2) = Eo = const. (45) 

which reflects the conservation of energy at the microscopic level. 

Remark 1: The evolution is described by Eqs. (41) and (42) together with Boltzmann equations for 
the fluid Eqs. (15) and (16). In the next section, we shall introduce the "average assumption" which 
will then enable us to ignore the equation for the fluids and to replace them by Eqs. (43) and (44). 

Remark 2: Recalling that A2{t = 0) = 0, it is straighforward using Gronwall lemma to check 

the consistency of the scaling hypothesis, i.e. A2 = 0{a.), from the above evolution equation (42). 
Indeed, F3 remains of order 0(1) (upper bound denoted sup(F3)) and Fi is strictly positive since 
Fi = would imply that: 
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which is precluded by the relations: 

> V{t) (1 - a) pj^^^iv, t) = p_/2V^, - + - V,).t) 

yv<V{t) {l-a)pt^,^{v,t)^pt^^j{2Vt~v + a{v-Vt),t) ^ > 

Therefore A2 is bounded by Qsup(F3)/inf(Fi) i.e. A2 = 0{a). 

Remark 3: In our first paper [10], we have considered the case a = 0, i.e. the thermodynamic limit 
for the piston, and discussed the evolution described by A2 = and: 

= F2{V) (48) 
= -MVF^{V) (49) 
= MVF+{V) (50) 





1 dV 




7 ~dt 


1 


d{E-) 


7 


dt 


1 


d{E+) 


7 


dt 



6 Adiabatic evolution for short time 

We want to investigate the evolution of the piston under the initial condition where X(t = 0) = Xq, 
V{t = Q) =0 and the fluids on both sides of the piston are in equilibrium, described by Maxwellian 
distribution of velocity, with temperatures T^. The initial conditions are such that |Tq" — Tq~| = 0{1) 
and ^ Pa \ = C'(l). Since a ^ 1, the perturbation approach shows that in a first stage, we can 
restrict our analysis to order zero in a and thus we recover the results of our previous paper (where 
a = 0), except that they will now be valid only for finite time. In this first stage, the motion of the 
piston is a deterministic (no velocity fluctuations, A2 = 0) adiabatic (no heat transfer between the 
compartments) evolution towards mechanical equilibrium. This first stage ends when p~ —p~^ = 0{a) 
(but T~ — T+ is still 0{1)). At this point, the term F2 which appears in the evolution of V, Eq. (41), 
becomes of order a and the first order terms of the perturbation approach must now be taken into 
account. Let us note that if we introduce the "average assumption" (see below, Assumption 2), the 
state at the end of this first stage is given by the final state derived in [10] for a = 0: 

Pid = Po + Oia) (51) 

■Oia) (52) 



-) 


1 To 


Xad 


N-) 


L 


N+ ) 


1 To 





where: 



= (ik) To {l-^)+0{a) (53) 



ksTo=i^) Po = ^, N^N-+N+, (54) 



N J N 



and Xad is the solution of 



'^)^L->/(^)(i-W^y|c (55) 



where the constant C is related to the initial conditions according to 



C = VT-(0)Xo - x/r+(0)(L - Xo) (56) 
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7 Slow relaxation towards thermal equilibrium 



The perturbation approach developped here (at order 1 in a) intends to reach the long-time behavior 
of the piston motion in the case where M ^ m. We claim that the fluctuations of the piston velocity 
V enter the scene only in a second stage, of time scale 0{l/a). once the pressure difference has become 
of order 0{a). In order to prove this assertion, we have to develop a boundary-layer-type perturbation 
approach[20]; indeed, at the end of the fast relaxation towards mechanical equilibrium, the standard 
perturbation approach (with time variable t) becomes singular and fails to give access to the further 
evolution of the system. The relevant time variable to be used in order to reach the second stage of 
the evolution is the rescaled variable: 

T = at (57) 

When the evolution is described in terms of this rescaled time, the first stage collapses into a boundary 
layer (r <C 1) whereas the focus now bears on the second (previously asymptotic) stage. We introduce 
a rescaled velocity V{t), describing a slow motion of the piston: 



V{t) = a V{t) 



with V{t) = 



dX 
dT 



(58) 



The motion is now driven by fluctuations and the velocity fluctuations play a crucial role; they allow 
to define the piston temperature T^: 



A2{t) = a A2(t) 



with A2 = 



2m 



(59) 



Finally, mechanical equilibrium has been reached in the first stage of the evolution and only fluctua- 
tions around mechanical equilibrium are to be observed in the second stage. We thus introduce: 



{p- -p+){t)=aIl{T) 



(60) 



In this second stage, the rescaled quantities are of order 1 and this stage ends when they become of 
order a. At this point, we should then consider the corresponding equation to order a^. We now 
investigate the consequences of the perturbation approach to order 1 in a. Equations (41)-(44) can 
be written in terms of the rescaled quantities: 



( adV^ n 
7 ar 2m 



(61) 



a dA2 



-4A2F1 + F3 



7 dr 

Another set of equations describes the evolution of the gas energies (at lowest order in a): 



1 d{E^) 
W dr 

1 d{E+) 



2m 



A 



V [F^- + A2F0-] + m ( [4A2i^f 



f; 



2m ( ^1 V [F+ + A2F+] - m 



(62) 



N+ dr "' \N+ J 

Consistency of the perturbation approach then requires to take the value of all the functions Fj at 
V = {as soon as V remains of order 0{a), i.e. V = C(l)), and therefore to replace Eqs. (62) by: 

1 d{E-) ^ A \ ~ _ (A 



N- dr 

1 d{E+ 
iV+ dr 



surf + m[—][4A2F,--Fi] 



2m 



A 



N+) 



^Pturf - 



m 



[4A2F+-F+] 



(63) 
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8 Slaving principle (consistency condition) 



Equations (61) show that V and A2 arc slaved to the slow relaxation of the gases towards thermal 
equilibrium appearing in the p^-dependence of the Fj: 



(64) 

~ ^ ^3(0) 

' 4Fi(0) 



We call consistency condition such a lower-order resolution of evolution equations of the general form 
a dA/dr = r.h.s.{T) — 0(1), leading to the instantaneous relation r.h.s.{T) = 0, up to terms of order 
0{a). Indeed, it merely follows from a term-wise identification of the expansion in powers of a. given 
that A actually is a non trivial function of r when a 0[18]. It can also be termed a quasi-steady- 
state approximation (as for instance in the context of enzymatic catalysis [19]). We shall make further 
encounter of this argument. As mentionned above, the consistency of the lower order perturbation 
approach also requires to set V = in functions Fj. 



9 Average assumption 

At this stage, to simplify our analysis, wc introduce the following: 

Assumption 2 (average assumption): The temperatures and the densities at the surface of piston co- 
incide at order 1 in a with the average energy and density in the fluids in the respectively left/right 

compartments : 

a) T^urf = T^ where iV± fcs T± = 2(£;±). 

N- , N+ 



b) Psurf=P wliere p = 
From Eq. ( 45), we get: 



p+ = 



A{L-X)' 



and then p^^^^ = p"^ = p^ksT^- 



or equivalently 



Therefore: 



N-ksT- + N+kBT+ = 2Eo - M(A2 + V^) 



Xp- + (L-X)p+ = ^-^(A2 + y^) 



p = po + a 



p+ = Po - a 



X 



M 



1-- n-— (A2 + ay 



AL' 



M 



where by definition pq = 2Eq/AL, and thus: 



dp^ 
Ik 



n+^(A2) + al/2) 



= 0{a) 



(65) 
(66) 

(67) 
(68) 

(69) 



10 



Moreover, from Eqs. (63) wo have: 



dr 



A 

TF 



Vp- + 2m ( [4A2Ff +F3-] + 0{a) 



dT+ A 



(70) 



Vp+ - 2m 



A 



[4A2F+ + F+] + 0{a) 



On the other hand, from Assumption 2b (average assumption), we have: 



N-keT- 

which together with Eqs. (70) and (71), yields: 



= -3p-V + 





1 


dp 






~d^ 




1 


dp+ 






Ik 



2m 



N+kBT+ 
A{L - X) 



[4A2Ff -F3-] + 0{a) 



2m 



3p+V - ^[AA2F+-F+] + Oia) 



From Eq. (64) follows that: 



4A2Fr - F3- = 4A2F+ -F^ = Y^ {F3F+ - F+F,-) 

On the other hand, from Eqs. (67) and (68), we have: 

p-ksT- = p+kBT+ + 0{a) 



(71) 

(72) 
(73) 

(74) 
(75) 



Therefore, since dp^/dr is of order a, the consistency condition of the perturbation approach yields 
from Eqs. (72), (73) and (74) the final relation: 



V = 

h _ 

2m ~ 

Ao = 



fAL\ fF,F+-F+F,- 



m 



+ 0{a) 



2m f AL 
~3 



(76) 



_F3_ 

4Fi 



+ 0{a) 



10 Dimensionless variables 

As usual in perturbation methods, it is convenient to introduce dimensionless variables. Prom Eq. 
(14) and Assumption 2a, wc have: 



Let us define: 



N-ksT- + N+kBT+ = 2Eo - M(A2 + V^) 
2Eo fAV 



ksTo = 



N 



N 



Po 



It is natural to introduce the dimensionless variable (, defined by: 

N^ksT^ ^lil±2^) iNkBTo- Ma{A2 + aV^) 



(77) 
(78) 

(79) 
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We thus have from Eqs. (67) and (68): 



N+kBT+ - N-knT- 



2^[ALpo - Ma{ A2 + aV^)] 
p+ A{L - X) ~ p- AX 



= PqAL- aAL 

-2poAX + aAX 



f2X \~ 2M ~ -2,' 



Therefore: 



2€ 



ALpo 



L J Po 



X 



L-1 AL 



M ~ ~ ' 

— n A2 



(80) 
(81) 

(82) 

(83) 

(84) 



The variable ^ remains C(l) in tlic second stage of the evolution, hence the consistency of the per- 
turbation approach requires to truncate the previous equation at lower order in a, which yields: 



According to this the consistency condition, we have: 

di__Y_ 

and from Eq. (76), we obtain our final equation: 



(85) 



(86) 



(87) 



11 Maxwellian identities 

Let us introduce at this point the following: 
Assumption 3 (Maxwellian identities): 

a) The relation between the functional Fi and F3 obtained when are Maxwellian remains valid 
here at order a, i.e. 



Fi'iV) = ^^^^ F^{V) - VF^{V) + 0{a) provided V = 0{a) 
b) = F^{V = 0) coincide at order a with the value given by Maxwellian distributions, i.e.: 



0{a) 



(88) 



(89) 



From Assumption 3a, we have: 



01. 

FiF+ - F+F,- = (T- - T+)FrF+ 



and thus: 



dT 



/ 2A \ F-F+ 



(90) 



(91) 
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kB {T-F^-T+F+) 



A2 = — 

^ = 'B(r--T+) 



En 



ALJ k%T-T+ 



(92) 
(93) 
(94) 



We remark that whatever are the distributions p^{v, t), then as long as Assumption 3a is satisfied, we 
have V > if and only if T+ > T~, i.e. the piston moves in the direction of the warmer side, which 
is typical of heat transfer since p+ = p~ + 0{a). 

To obtain an explicit equation for the evolution, we shall now use the Assumption 3b above, 
together with the fact that p'^kBT'^ = Po + 0{a). We thus have: 



F^F+ 



AL 



VkE Vt^ + Vt+ 



and Eq. (87) leads to the main result of our paper: 



dT 



1 


1 




FY' 




3 


L ^ 


/ N ^ 


V rmr 





N 



N+ 



N 



N- 



Finally introducing the dimensionless variable s by: 



- ...^ . 2N- 2N+ 



N 



N ^0 



(95) 



(96) 



(97) 



we have: 

Main result: under Assumptions 1, 2 and 3, the evolution of the system in the second stage is described 
in terms of the dimensionless variables by 



ds 



N 
2N+ 



(1 + 20 



TV 
2N- 



(1-20 



(98) 



together with the "initial condition" (in fact the matching condition of the boundary-layer-type per- 
turbation approach) which is from Eq. (85): 



as = 0) = 



1 x„ 



(99) 



with Xad being the position of the piston at the end of the first stage, Eq. (55), i.e. at mechanical 
equilibrium. Integrating Eq. (96) with Eq. (99) yields the evolution ^ = C('S) from which we obtain 
the position and the temperatures: 



X 



2N± 



(100) 
(101) 
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together with V, A2 and 11 obtained from Eqs. (92)-(94) and Eq. (95), i.e.: 



3V 7rm V 



A2 = -2-VT-T+ i.e. = VT+T- 

2m 

Eo\ fT+-T- 



n 



(102) 



,.--1 



Prom Eq. (96), the second stage of the evolution described here will proceed until: 

N+ -N- 



which implies: 



^=0 



2N 



0{a) 



NkB 



f 



Eq 

~" 2m mN 

V = 0{a) 



+ 0{a) i.e. Tf =Tf +0{a) 



(103) 

(104) 
(105) 

(106) 
(107) 



n = 0{a) 



(108) 



At this point, corrections in have to be included; the relevant perturbation approach to get the 
further evolution would use the rescaled time T2 = o?t and consider also A3 and A4, and so on for 
the successive stages of the relaxation towards complete equilibrium (Maxwellian distributions for the 
gas particles). 



12 Remarks 

1) Our main result Eq. (98) gives an equation totally independent of the parameters of the problem, 
except for the ratio 2N^ /N , which will be 1 in our simulations. 

2) The expression (102) for V and the equation (96) describing the evolution, coincide with the 

equations obtained in [9] using qualitative arguments. 

3) The equation (98) for ^(s) was also derived by J.L. Lebowitz [21] using either qualitative arguments 
or the fact that the evolution of the piston will be described by an Ornstein-Uhlenbeck process. 

4) The consistency between the short time solution (obtained in [10] with a = 0) and the solution 
for large time, here obtained by considering evolution equations with respect to rescaled variable r 
and restricted at first order in a, amounts to match the infinite-time limit t ^ 00 of the zero-order 
solution with the initial condition in r = of the rescaled perturbation large-time solution, which 
makes sense only for large R. 

5) Prom Eq. (102), we see that 11 > if and only if T+ > T~ . As we have seen above, it is precisely 
the condition for V > Q. Therefore work is delivered to the warmer side but more heat is extracted 
from this warmer side so that its energy decreases. 
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6) We should compare the expression obtained for the velocity Eq. (102): 

with the expression obtained in [7] for the stationnary state of the piston in an infinite cylinder under 
the condition that the pressures on both sides are equal: 

The velocity V is larger than V^t by a factor 16 /Sw. It is related to the fact that in the present case, 
the pressures are not equal but: 

p--p+ = ^ + Oia^)>0 (111) 



13 Numerical simulations 

In order to check the assumptions 1, 2 and 3 and to test our main results, in particular Eqs. (96)-(108), 
we have conducted a large number of simulations. We have considered a one-dimensional system of 
fixed length L. The mass of the fluid particles is m = 1 and the mass M of the piston is varied in the 
different simulations. Initially the piston is at the position Xq with velocity Vq = and a configuration 
of N~ and particles is taken at random from Maxwellian distributions with temperatures Tq and 
Tq". The system is then led to evolve according to the law of purely elastic collisions Eq. (1). The 
position X of the piston is recorded, together with the temperatures and pressures in the left and 
right compartments defined by the average (kinetic) energy and the equation of state ~ p^ksT^- 
To conduct these numerical simulations, one is confronted with the following dilemma. In order 
to reduce the fluctuations to a minimum to obtain significant results, one should take very large. 
This implies that M should also be very large in order to have a strong damping. However since the 
second stage proceeds on a time scale r = 2t/M, this implies that one has to follow an enormous 
number of collisions, which needs a very long computer time. For this reason, to study the first stage 
of the evolution, we have considered up to 2 millions particles, but we were forced to reduce drastically 
this number of particles to investigate the second stage, i.e. to approach thermal equilibrium. In all 
simulations presented below wo have taken: 

N 

kB = l, m = l, L = 60, N-=N+=N= — 

Xo = 10, Vo = 0, Tq- = 1, T+ = 10 i.e. = 5p+ and Po = ^ 

For these initial conditions, the end of the first stage, i.e. mechanical equilibrium, is characterized by 
Eqs. (51)-(55) which give: 

X„rf = 8.42, T-, = 1.54, T+=9.46, p" = p+ = 0.1833 TV (112) 

Wo conchidcd in [10] that in the first stage, the motion is deterministic and depends strongly on 
R = N/M if i? < 1, but become independent of if i? > 4. In particidar for i? < 1, it is weakly 
damped and the period of oscillations we have computed coincide with the period obtained from 
simulations and from hydrodynamics assuming adiabatic oscillations. On the other hand for R > 10, 
the motion is strongly damped and is independent of and M as soon as M is large enough. 
Moreover at the end of the first stage, after mechanical equilibrium has been reached, simulations 
shows that oscillations start to appear. These oscillations were interpreted as a consequence of the 
fact that the velocity distributions in the fluid were not Maxwellian. 
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In Figure 1, we have plotted the first part of the evolution, i.e. the adiabatic approach towards 
mechanical equilibrium, in function of time t for = 2.10^, M = 10^ (i.e. R = 0.2, weak damping) 
and for A^''^ = 3.10'', M = 3.10^ (i.e. R = 10, strong damping). From the simulations the values 
obtained at mechanical equilibrium are: 

Xobs = 8.33 ±0.05, T-^^ = 1.52 ±0.04, = 9.48 ± 0.04, p"^ = p+ = 0.183 iV (113) 

in good agreement with the analytical values Eq. (112). Other graphs for this adiabatic evolution can 
be found in [10]. 

In Figures 2-8, we have plotted the second part of the evolution, i.e. the approach towards 
thermal equilibrium in terms of the scaled variable r = 2t/M, obtained from our equations and from 
simulations with A'^='= = 3.10'* and M = 100, 200, 1000. The position in function of t is presented 
in Figure 2, where no significant difference can be seen for the different values of M. Figure 3 shows 
a zoom around t = and r = 16 for M — 200 and M = 1000, expressed in terms of the real time 
t = tM/2. It shows that on the time interval {t, t + 300) considered in this zoom, the evolution (fast 
local relaxation slaved to the slow evolution of the system) will tend to be independent of M for M 
sufficiently large (fluctuations around mechanical equilibrium) while on the large scale the evolution is 
scaled by 1/M, i.e. appears independent of M when expressed in function of r = 2t/M. In Figure 4, 
we present the temperatures together with the surface temperatures T"^^^ computed from the 
incoming particles at distance 1 from the piston. In Figure 5 the temperature of the piston is plotted 
using either Eq. (102) i.e. = VT+T' or the definition (77) with ksT^ = MA2; one should note 
that using the definition (77) implies a factor which is in our simulations 30000(T+ + T~ — 11) and 
thus the fluctuations are very large. The other scaled quantities of order 1, i.e. 11 = M{p~ — p+)/2 
and V = MV/2 (recall that to = 1) obtained from the simulations and our Eqs. (102) are represented 
in Figures 6 and 7. Finally, since the distributions of velocities of the gases do not remain Maxwellian, 
we compare in Figure 8 the previous evolution with the evolution starting from the (adiabatic) initial 
conditions (113) with Maxwellian distributions. We observe no significant difference. 



14 Conclusions 

Using a two-time-scale perturbation approach, we have shown that the evolution of a piston with large 
but finite mass proceeds in two stages. For R > 10, there is strong damping and in the first stage, at 
times t = C(l), the evolution is similar to the one previously discussed in the thermodynamic limit 
[10]: it is a deterministic, adiabatic evolution towards mechanical equilibrium. At the end of this 
first stage, adiabatic oscillations appear which are interpreted as associated with the fact that the 
velocity distributions of the gas particles are not Maxwellian. In particular, simulations show that in 
the adiabatic evolution: 

J-^.-1.8l^ (114) 

In the second stage, at times t = 0(M), a fluctuation-driven regime develops, leading to a relaxation 
towards thermal equilibrium. In this regime, the motion of the piston is slaved (i.e. it adapts on a 
time scale 0(1)) to the slowly relaxing (on a time scale 0(M)) asymmetry of the thermal fluctuations 
on each side. Moreover, since the motion of the piston is now stochastic, one has to introduce the 
temperature of the piston. Therefore the fluctuations in the piston motion produce a heat transfer 
from the warmer to the colder side which is larger than the work produced by the piston motion. In 
particular, in this stage, simulations as well as equations show that now: 

J_dr2^ ]_dx 

T- dt ^ X dt ^ ' 

In this regime, the piston is no longer adiabatic but the time involved to reach thermal equilibrium is 
of order M. 
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In the first stage, simulations are in very good agreement with the predicted values for the position 
of the piston and the temperatures of the gases at mechanical equilibrium. Moreover for weak damping 
(i? < 1) the observed values for the period coincide within a few per cent with the predicted values 
from our equations, as well as with the values obtained from thermodynamics assuming "adiabatic 
oscillations". However the observed values for the damping coefficient disagree strongly with the 
computed values. Unravelling the origin of the friction, wc must conclude that there are two different 
friction coefficients. A first coefficient describes the stabilization of the piston velocity in the infinite 
length situation, or for the case of finite length (and i? > 10) until the first recollision on the piston 
takes place. It is related to the shock waves propagating uniformly on both sides of the surface of 
the piston. The second coefficient describing the damping of the oscillations is associated with the 
rebouncing of these shock waves on the surface of the cylinder and of the piston. It appears to be 
related to the relative velocity of the shock waves with respect to the piston velocity [11]. 

For the second stage, as can be seen on Figures 2 to 7, our main conclusions are all verified by 
numerical simulations with high degree of accuracy. One should however observe that the agreement 
would be "perfect" if we would consider, instead of Eq. (98), the equation: 



The origin of the discrepancy could be traced to our Assumption 3b where we use the Maxwellian 
values for F^iV = 0). It is known [9] that at this stage the distribution functions f^iv) differ 
considerably from the Maxwellian distributions around v = 0. Of course this discrepancy could also 
come from the average assumption (Assumption 2); however from Figure 4, we see that the surface 
temperatures oscillate around the bulk temperatures (on a time scale t) and thus we expect this 
technical assumption to have no important consequences. 

Let us remark that if we compute the time auto-correlation function C{t) of the position X of 
the piston and the associated power spectrum S{u)) (the Fourier transform of C{t)), then from the 
simulation wc could expect a three-peak structure as observed in the equilibrium state by White et 
al. [15] recovering a standard result in hydrodynamics [22]. The central peak describes (according to 
the fluctuation-dissipation theorem) both 

the (slow) relaxation towards thermal equilibrium, described by a term ^ p-*/'^th jjj ^j^g tji^jg 
auto-correlation function C{t), with Tf?j ~ 1/a ~ M; the width of the peak is 0{1/M). 

- and the thermal fluctuations of the piston as described by A2(oo). 

The two lateral peaks, slowly evolving in time with X(t) and T^{t) correspond to the fast adiabatic 
and deterministic relaxation towards the instantaneous mechanical equilibrium, slaved to the values 
of X{t) and T='=(t). It also gives (fluctuation-dissipation theorem in this quasi-equilibrium state) a 
contribution to the fluctuations of the piston motion. The width of these peaks is independent of M. 

The perturbation approach can be carried on at higher orders; it would give access to the following 

stages of the relaxation, at increasing time scales 0(l/a^) ^ 0{M"), involving higher moments A2„_i 
and A2„. Maxwellian distributions for the velocity distributions inside the gases would be reached 
only in the final stage s of this relaxation (namely at times t ~ 0(M") with n large ). However 
the factorization property will be violated at these orders and we can not say what would be the 
consequences of this violation on the results obtained within our Boltzmann's equation approach. 
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Figure captions 



Figure 1: Adiabatic stage of the evolution for Tq = 1 and Tq = 10. 

a) strong damping: iV± = 3.10^ M = 3.10^, R = 10. 

b) weak damping: = 2.10^, M = 10^ R = 0.2. 

Figure 2: Approach to thermal equilibrium for = 3.10^ and M = 100,200, 1000 compared with 
the solution of Eq. (96) (the real time ist = tM/2). 

Figure 3: Zoom on the evolution of Figure 2 in function of the real time t for M = 200 (light curve) 
and M = 1000 (bold curve) around r = and r = 16. 

Figure 4: 

a) Evolution of the temperatures T±(t) for iV± = 3.10^ and M = 200 compared with Eqs. (96) 

and (101); 

b) Surface temperatures T^^j:{t) compared with bulk temperatures T^(r). 

Figure 5: Temper ature o f the piston {N^ = 3.10", M = 200) 

a) from = y/ T+T- (Eq. (102) with r± obtained from Eq. (96) (light curve); 

b) from = y/T+T~ with obtained from the simulations (dark curve); 

c) from the definition Eq. (77) with ksT^ = MA2 (stochastic curve). 

Figure 6: Scaled pressure difference Att = A[p~ — p+)M/2 computed from Eq. (102) with T"^ 
obtained from Eq. (96) and from simulations (with = 3.10", M = 200). 

Figure 7: Scaled velocities {N^ = 3.10", M = 200) 

a) from Eq. (96) (light curve); 

b) from Eq. (102) with T"^ obtained from the simulations (dark curves); 

c) average velocity obtained from the simulations. 

Figure 8: Evolution for = 3.10", M = 200 compared with the evolution starting from initial 
conditions Eq. (113) and Maxwellian velocity distributions. 
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